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We develop a simple thermodynamic model to describe the heat transfer mechanisms and gener- 
ation of acoustic waves in photoacoustic Raman spectroscopy by small particulate suspensions in a 
gas. Using Langevin methods to describe the thermal noise we study the signal and noise properties, 
and from the noise equivalent power we determine the minimum number density of the suspended 
particles that can be detected. We find that for some relevant cases, as few as 100 particles per 
cubic meter can be detected. 
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I. INTRODUCTION 

Optoacoustics as a field has venerable roots dating 
back to 1880 when Alexander Graham BellQ first found 
that periodic illumination of colored substances gener- 
ates sound. The process involves absorption of energy 
from a modulated light source by the atoms or molecules 
in the medium, producing a variable heat source which 
serves as a source for acoustic waves; the sound can be 
detected by microphones or piezoelectric transducers. 

Photoacoustic spectroscopy, as it is also known, has 
become a valuable spectroscopic tool with extensive use 
in the study of absorption spectra and electronic spec- 
tra of substances in both condensed and gaseous phases 
H H H 0. The type of atomic or molecular modes be- 
ing studied determines the frequency of the incident light 
used, which can be in the visible range for electronic ex- 
citations or in the infra-red for vibrational modes. 

Photoacoustic techniques were first applied to Raman 
spectroscopy in 1979 by Barrett and Berry 6], to cre- 
ate the nonlinear spectroscopic technique they called 
photoacoustic Raman spectroscopy (PARS). It was sub- 
sequently applied to liquids and to gaseous trace 
analysis [g. Resolution much higher than the natural 
linewidth of Raman transitions has been achieved using 
narrow linewidth pulsed lasers 0- 

Numerous researchers in the last couple of decades 
have routinely applied photoacoustic methods to deter- 
mine structural and thermal properties in solid state 
physics 10, [Till, t o conduct non-intrusive studies of 
nanostructures |l2l |. to analyze internal structure and 
dynamics of molecules in chemistry |l3j . and as a non- 
invasive way of characterizing wavelen gth dependent op- 
tical properties of biological samples [Tjl- But almost 
all such applications have involved linear infrared exci- 
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tations, and very often the photoacoustic measurements 
were a supplement to separate Raman spectroscopic anal- 
ysis. In contrast there has been limited literature on the 
usage of PARS as a tool by itself for spectral studies. 

However as we will show here, PARS can become a 
very effective and sensitive means of spectral detection 
of suspended solid contaminants in air. The immediate 
motivation for our study arises from the urgent neces- 
sity for reliable methods to detect trace amounts of mi- 
crobial spores like anthrax in the current context of de- 
fense against bio-terrorism. Methods based on coherent 
Raman spectroscopy are already being actively pursued 
|15| where the spectroscopic signature of a characteristic 
molecular constituent is sought. An optoacoustic method 
based on Raman spectroscopy can take advantage of the 
same signature but offer an alternative and experimen- 
tally simpler detection scheme. In fact the first appli- 
cation of infra-red photoacoustic spectroscopy for iden- 
tifying bacterial spores was reported just a few months 
ago jig where the constituents in a mixture of diverse 
specimen were identified correctly with 100% accuracy. 
Applications are viable also in a broader context since 
acoustic detection of stimulated Raman emissions using 
the powerful lasers available these days can lead to more 
sensitive trace analysis. 

It is worth noting that application of photoacoustics 
to highly sensitive detection of trace amounts of pollu- 
tants in gaseous medium was originally responsible for 
its tremendous growth and revival in the 1970's. Tradi- 
tionally such applications have almost always relied on 
infrared single photon absorption by the pollutants. But 
using Raman scattering provides several advantages: (i) 
the incident light can be in the visible range in which both 
air and water are essentially transparent while they ab- 
sorb in the infrared; this is particularly relevant for bac- 
terial spores because most of the non-signature molecules 
are water, (ii) adjusting the difference of the two input 
fields allows tuning to different vibrational modes not 
available for direct infrared excitation, and (iii) in stimu- 
lated Raman scattering the signal depends on the prod- 
uct of two laser intensities, and so it can be maximized 



using strong lasers. 

Moreover new innovations in Raman spectroscopy, 
such as STImulated Raman Adiabatic Passage (STI- 
RAP) [3 or using chirped pulses allows complete pop- 
ulation transfer from the ground state to the first molec- 
ular vibrational state. Subsequent relaxation of the ex- 
cited state can generate sound waves. In fact experimen- 
tal evidence already exists , that a coherently excited 
molecular gas generates intense sound waves when the 
incident electromagnetic fields are close to Raman two- 
photon resonance. 

In this paper we will consider a model for a PARS ex- 
periment to detect small (~ 1/im) suspended solid par- 
ticle in air. We determine the signal strength by consid- 
ering the basic thermodynamic processes involved and 
describe the acoustic disturbance by a wave equation. 
To describe the noise in the system we use a Langevin 
method which has the advantage of being able to include 
possible additional sources of noise with little extra effort. 
From the expressions for the signal and the noise we de- 
termine the noise equivalent power (NEP) and from that 
the minimum detectable particle density. We describe 
our model in Sec. II, and then in Sees. Ill and IV we 
consider the thermodynamics of the energy transfer from 
the fields to the photoacoustic cell. In the following two 
sections V and VI we describe the acoustic waves and 
evaluate the thermal noise and the NEP using second or- 
der Langevin equations, and then in Sec. VII we bring 
all the elements together to determine the minimum de- 
tectable number density for the specific example of an 
anthrax spore. 



II. THE MODEL 

The configuration we consider, shown in Fig.^ is typ- 
ical of most photacoustic experiments where the driving 
lasers pass through a gas/air-filled cylinder which may 
contain the Raman active medium of interest in the form 
of suspended particles. The detector is located at a po- 
sition orthogonal to the direction of propagation of the 
laser beams. Mirrors may be placed at the ends of the 
gas cell for multiple passes of the lasers. 

The Raman scattering needs to be Stokes- type because 
energy needs to be deposited in the medium, and neces- 
sarily be stimulated in nature because of the stronger 
signals produced. Thus there are two input beams: the 
pump beam with energy hi^p and a Stokes beam with en- 
ergy < hvp. The Raman scattering scheme is shown 
in Fig. 121 with a and b being the upper and the lower 
levels and energy difference is tioja — fiLUi, — hujab ■ Raman 
resonance is achieved when i^p — ly^ — i-^ab- 

Each of the spores we wish to detect has volume Vs 
and contains Ns molecules; out of these a fraction fj^ 
are the Raman active molecules we call R, and the re- 
maining 1 — /ij we take to be predominantly composed of 
another species of molecule W, usually water for a bac- 
terial spore, that is effectively non-absorbing in either of 
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FIG. 1: Schematic configuration of our model for optoacoustic 
detection of suspended particles (indicated by small circles) in 
air contained inside a photoacoustic cell (shown as a cylinder) . 
The drive and the Stokes photons carry energies hvp and hva . 



the input fields. The Raman scattering deposits energy 
in the gas cell whereby heat is generated that serves as 
the driving force for acoustic waves. The heat-energy 
transfer happens in two steps: 

Step 1. The excitation energy Tiujab of the molecules R 
will be converted to heat energy via collisions R ^ W, 
R ^ R and W ^ W among the molecules in a particle. 
This will depend on the Raman scattering cross-section 
and the radiative and coUisional damping that happens 
within a spore. 

Step 2. The energy gained by molecules in the spore 
will be transferred to the gas molecules in the tube. Since 
the spore is much larger than the gas molecules, this can 
be modelled by kinetic theory. 

If we define Hn to be the energy absorbed due to Ra- 
man scattering per unit volume per unit time, then the 
energy absorbed by a spore per unit time is HrVs- A 
fraction rj of this energy is converted to translational en- 
ergy or heat of the surrounding gas molecules within the 
time scales of acoustic modes, and therefore that is the 
fraction available for the acoustic signal. The spores will 
be randomly distributed in a gas cell, and we define ps to 
be their average number density. The amount of energy 
available for acoustic signals per unit volume of air per 
unit time is then given by: 

H = PS X T] X HrX Vs (1) 

In order to determine the minimum detectable ps, we 
will have to find the minimum energy density H that can 
produce a detectable signal, and for this we first need to 
estimate all the other quantities in the above expression 
as well as the noise level. 



III. RAMAN INDUCED HEAT ENERGY 

In the Raman scattering process, the pump beam loses 
photons by exciting the molecules and their subsequent 
de-excitation produces Stokes photons, and the energy 
difference is deposited in the medium. For each photon 
converted the deposited energy is hvp — TiVg. The total 
energy deposited is determined by the number of photons 
gained by the Stokes beam, which has an exponential 



3 



dependence on the interaction length z, 



n,(z) =r^,(0)e»^ 
The gain coefficient g is given by ^ 



(2) 
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where Vg — i x 10*ms^^/n(i^s) is the velocity of the 
Stokes photons with n{vs) being the refractive index in 
the medium; Sv is the hnewidth for the Raman transi- 
tion; N is the density of the active medium. The last 
factor is the fractional population difference of the upper 
and lower atomic levels for thermal occupation proba- 
bilities. The gain is relatively small over the interaction 
lengths L which is on the scale of a spore-size for trace 
occurrences, so we may write e^"^ ~ 1 -|- gL. This also 
means that the fractional change in the input beam in- 
tensities is small so that we can treat the average values 
of Ip and Is as constants during the time of interaction. 
Thus the change in the intensity of the Stokes signal is 
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(4) 



where we defined the intensity-independent gain factor 
G = g/Ip- By definition, intensity = power/area, there- 
fore dividing through by L we have a relation for the rate 
of change of energy of the Stokes beam per unit volume. 
Each photon converted deposits energy hvp — TiVs in the 
medium, so the rate of energy absorbtion per unit volume 
of the medium is given by 
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(5) 



There are three possible fates for this energy: (i) it can 
be lost by radiation, (ii) energy from one molecule can 
be transferred to excite another molecule via collisions 
or (iii) the energy will be converted to translational en- 
ergy via collisions. The last process is the one which will 
directly lead to heating of the system. Since the exci- 
tations are vibrational in nature the second process will 
also cause heating as well after a few collisions. We de- 
note the radiative decay rate by which corresponds to 
the first process of loss, and we denote the collisional loss 
by Fc which corresponds to the last two processes. With 
these considerations we can now relate the heat absorbed 
per unit volume of the medium to the Raman gain factor 
and the input laser intensities 
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IV. HEAT TRANSFER EFFICIENCY 



We next need to determine the part of this heat that 
is actually available as acoustic energy. This is parame- 
terized by T] which depends on how much of the Raman 
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FIG. 2: Stokes-Raman scattering. The driving field has fre- 
quency Vp and the Stokes field Vs , and they induce transitions 
between levels a and h via an intermediate state c which could 
be virtual but may coincide with a real energy level. 



energy absorbed by a spore can be transferred to the sur- 
rounding gas in the optoacoustic cell on the time scale of 
acoustic modes. As we noted in the previous section the 
heat is lost from the spore through radiation and through 
collisions with gas molecules. 

We first calculate the collisional energy transfer rate. 
Assuming a spherical particle of radius Ts , the rate of col- 
lisions Nc of gas molecules with the surface of the particle 
is determined by 



Pq X AttK 



Nr. 



X 2mgUg 



(7) 



with mean velocity 
m 



Ug of the gas molecules given by 
~ ^iiTigkTo; at standard temperature and pres- 
sure (STP), Nc 4: X 10^^ collisions/second. The tem- 
perature increase of the spore due to Raman absorption 
is given by 



c^AT = fRNah{Vp - Vs) 



(8) 



where Na is the Avogadro number and the molar 
specific heat. Most of the molecules in a bacterial 
spore are usually water (for example in anthrax spores 
fn ~ 0.1) and so we use the specific heat of water 
Cy = 4.2 J/mole/°K. Then for typical vibrational fre- 
quencies of ^ 10^'* Hz and Raman active fraction in the 
range = 0.1 — 0.5, the spore will gain a temperature of 
100° — 500°K. Thus the increase in temperature Ts of the 
spore after Raman scattering is of the order of the initial 
temperature of the gas (~ 300° K), and since the water 
and the diatomic gas molecules have similar masses, the 
energy transfer in each collision will also be of the same 
order ^ kTo ^ kTs- Therefore the rate of heat transfer 
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from the particle to the gas is given by 

c„r, ~ -NckTs X ^ = -^RTs, (9) 

and the time scale of the collisional heat loss is r ~ 
CvNs / RNc- The number of water molecules in a micron- 
size particle is Vs x iVo ~ 10^^. Thus t ^ 10"^s < the 
typical acoustic modulation times of ~ 10^^ s. 
The rate of heat loss by radiation is given by 

UB X X AtttI W (10) 

where ctb = 5.670 x 10"® Wm"2j^-4 jg ^-j^g Stefan- 
Boltzmann constant. At the upper limiting temperature 
of 1000°K for a micron-size particle the radiative flux is 
~ 10~^W. The energy deposited in the spore is given by 
Nsh{vp — Vg) ~ 10~® J, so the time scale of radiative loss 
is of the order of 0.1 s, which is much slower than the 
collisional loss rate. Therefore we can assume that on 
the time scale of the acoustic waves, essentially all the 
heat from the spores is transferred to the gas through 
collisions, and so we can take -q — 1. 



The sound velocity in the medium is = Pqj/po with 
Pa — *■ equilibrium density, V — > the volume of the gas cell 
and 7 — Cp/Cv — > the ratio of specific heat at constant 
pressure to that at constant volume. The macroscopic 
damping rate of the acoustic modes is denoted by T and 
the external driving force is F{f,t). We have added a 
random force R(t) to allow a Langevin treatment of the 
fluctuations. The calculation of the noise in terms of the 
displacement variable has the intuitive appeal that the 
driving terms carry dimensions of force. 

The process of heat transfer on the other hand is best 
described in terms of a thermodynamic variable and we 
will write an equivalent wave equation for the pressure 
variation p which is related to Eq. (|11|) by 

Newton's law Pottt = — Vp and p = — poc^V^. (12) 

Before we write the equation we determine the form of 
the driving force. For fixed total volume of the gas in the 
cell, basic thermodynamics relates the change in heat to 
the pressure differential 



V. ACOUSTIC SIGNALS 

The absorbed heat will cause a pressure wave which 
will radiate out from each spore. Considering that the 
time for a laser pulse to traverse the optoacoustic device 
is very short compared to the time scale of sound prop- 
agation, the heat absorption and subsequent release is 
simultaneous for all the spores. So the wavefronts from 
each spore will be in sync and will add constructively to 
cause a macroscopic expansion/contraction of gas in the 
optoacoustic device. This allows us to describe the acous- 
tic waves in the tube by a macroscopic wave equation for 
which H serves as the source. 

When the density of the spores is sufficiently high to 
give a statistically uniform distribution, the wavefronts 
from all of them will merge to yield acoustic waves which 
bear the symmetry of the beam profile, typically cylin- 
drical; on the other hand if the density is such that the 
number of spores in the tube is of order unity, the waves 
will bear the symmetry of the spore which is roughly 
spherical. For an estimate of signal and noise we need 
not delve into the specific nature of the acoustic modes, 
except to note that only the lowest modes are signifi- 
cantly excited. 

There is no net transport of matter for weak acous- 
tic modes, so the waves can be described by relative 
displacement of infinitesimal volume segments in the 
medium; we denote the displacement from equilibrium 
of the air at distance r from the source by ^(r, t) and 
the associated variation from the equilibrium pressure by 
p(r, t) = P{r, t) — Po- The propagation of the disturbance 
is then described by 



dQ = dE^ —fT^dP ^ dp^ ^'^ ' dQ (13) 
rt V 

Earlier we defined H as the heat generated per second per 
unit volume of the gas due to the incident light so that 
dH = ^dQ/dt, it carries dimensions of power/ volume 
— pressure/time. This heat drives the acoustic waves 
and we can now write the wave equation for the pressure 
variations 

Here the random term Rp{t) is not a force, but it can 
be derived from the force R{t) in the displacement equa- 
tion. The first step to the solution of the wave equation 
is to eliminate the spatial derivative by doing an expan- 
sion in terms of the normal modes of the corresponding 
homogenous equation 

p(r , t) = Y, A, {t)p, (r) ; (^V^ + p, (r) = (15) 

where the mode functions Pj {r) and the mode frequencies 
LOj are determined by the symmetry of the system and the 
boundary condition of vanishing velocity at the cell-walls. 
The displacement ^{r,t) has an identical expansion with 
the same frequencies but different mode functions ^jir) 
and amplitudes qj(t). The modes are orthogonal and 
normalized to 

J dVp,{r)p,{f)^VS,j. (16) 
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VI. SIGNAL, NOISE AND NEP 



to get the noise spectrum of pressure-wave amplitude 



The signal strength is determined by the mode ampli- 
tudes of the pressure variations that satisfy 



The signal damping rate Tg , which determines the profile 
of the acoustic signal, is determined primarily by macro- 
scopic thermal and viscous losses. After doing a Fourier 
transform, Aj{t) = J duje~'^'^*Aj{uj) etc., we obtain an 
expression for the mode amplitudes in the frequency do- 
main 



Aj{uj) 



iLu{j-l) J dVpj{f)H{f,uj) 



(18) 



To calculate the noise we use the normal mode expan- 
sion of the displacement in Eq. and write the resul- 
tant second order Langevin wave equation for the mode 
amplitudes qi{t) as two coupled linear equations 



du 



dt 



^^[F'it)+Rm ^ = (19) 



with F'(t) = -uj]qj{t) -I- ^ / dV^j{r)F{r,t). The damp- 
ing rate r„ of the noise amplitudes is mainly determined 
by the characteristics of the detector and should be dis- 
tinguished from the signal damping. A time-frequency 
Fourier transform of the first of the two equations, along 
with the property of the noise that {R{t)) = 0, deter- 
mines the power spectrum of the velocity distribution 



{Uj{uj)Uj{uj')) = 



(20) 



Assuming white noise, the noise is delta-correlated 
{R{uj)R{u}')) cx 5(u} — uj'). The Weiner-Khintchine theo- 
rem gives the correlations for the noise force in the time 
domain, which define the diffusion coefficient D, 



{R{t)R{t')) ^ 2D5{t - t') 



(21) 



An inverse Fourier transform and the equipartition theo- 
rem for each mode in thermal equilibrium, pQV{u'j{t)) — 
kT, yields the diffusion coefficient 



{Ui{t)ui{t')) 



D 



-r„(t-t') 



D = paVTnkp2) 



Now we take a Fourier transform of the second order 
equation in the displacement mode amplitudes 



qj{uj) 



F'{uj) + R{uj) 



(23) 



and use expression (|22|) for the diffusion coefficient and 
the relation between and p in Eq. 112|) 



'ipoCLOjqj 



(24) 



poc^u^TnkT 



(25) 



The noise equivalent power (NEP) is defined to be 
the input power necessary to produce a signal amplitude 
equal to the noise amplitude [23 ■ In order to estimate the 
NEP we note that the signal is typically generated in the 
lowest mode loq — Q corresponding to an uniform excita- 
tion of the gas. The noise is also dominated by the lowest 
modes for weak excitations. For the very lowest mode the 
noise amplitude vanishes, so we take the amplitude of the 
next higher one toi as a measure of the magnitude of the 
noise. Thus the noise equivalent power is approximately 
given by the value (indicated by subscript NEP) of the 
input power VH for which |^o('^)P — |^in(^^)P- Using 
the expressions from Eqs. (fH^ and (^5)1 we get 



\VHNEp{uj)f 



where we also used the fact that acoustic modulation 
frequency is typically much smaller than the mode fre- 
quencies to <^ OJj^Q. 



VII. MINIMUM DETECTABLE DENSITY 

The minimum detectable number density of spores 
Ps can be estimated by setting the input power den- 
sity H = HpfEQ, the noise equivalent power density in 
Eq. 1^. For the purpose of numerical estimation we will 
use the specific example of an anthrax spore for which the 
Raman active molecule is dipcolinic acid (DPA) which 
constitutes 17% of the weight, the rest being mainly wa- 
ter. Since some of the properties of DPA are not easily 
available, for those we use the values for benzene an or- 
ganic molecule similar in structure. 

The dimensions of gas cell are optimized for the maxi- 
mum interaction with the input beams, for a typical cell 
length of ^ = 0.1 m, the optimum radius for focussing the 
beams into the cell is given by r = ^/Xl/^^ ~ 1.4 x lO"'' 
if we use a wavelength in the visible range A ~ 7 x 10^ 



so that the cell volume V ~ nr'^l ~ 10 



TO 



We will 



assume that by placing multiple microphones all of the 
generated signal can be detected, but this can be easily 
relaxed by including a geometrical factor which involves 
the ratio of detector surface area to the cell surface area. 

At standard temperature and pressure, Tq = 300° K 
and Pq = 1 atm. and po — 1.3kg m^'^; for air 7 = 1.4 be- 
ing mainly diatomic gases. We take the dominant noise 
mode to be at the resonant frequency of the detector, typ- 
ically a condenser microphone, uji ~ cj™, and a reason- 
able value pol| is 4 X lO'^Hz, the damping of the mode has 
the same order of ma gnitude r„ = 5 x lO'^Hz assuming 
a quality factor = u/T = 0.8. The optimum modu- 
lation frequency is determined by the thermal damping 
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time w ~ Fs ~ 100s ^. Using these values we find 
Hnep{oj) ~ 5 X 10"^Wm-3si/2 

In order to calculate the Raman gain factor G, we 
use optical frequencies for the input beams i'p[s) ^ 
27r X 4 X lO^^'s"^. The number density of DPA molecules 
in an anthrax spore is 4 x 10^^ molecules m^"^. We 
use the Raman cross-section for Benzene ^ — 32.5 x 
10~'^'*m^/sr/molecule. A typical linewidth of a Raman 
line for Benzene is about 5v ^ 6.45 x lO^^Hz. Using 
these values we get for the gain factor 

G = 8.5 X lO^^VW^^ (26) 

Finally noting that an anthrax spore has dimensions 
1x2x1 ^TD? we can use Eq. to find the minimum 
detectable number density of anthrax spores in air 

HneP^/Ts 3 X 1025 _3 

PS ^ 77 71- = ™ 

7JX HrX Vs Ipis 

where we multiplied by the square root of the width of 
the acoustic signal since the NEP is not really "power" 
but power / \/ frequency . With powerful lasers available 
today intensities of /pjs) ~ lO^^Wm^^ are reasonable, in 
which case one would get ps ~ 30m~'^. 

Note that so one cannot indefinitely increase the laser 
intensities because at intensities of lO^^Wm"^ there is 
cascade breakdown of air at STP. While cooling the 
gas cell seems to be easy way to reduce thermal noise, 
due to condensation effects temperatures cannot be low- 
ered too much. Also it is much simpler and realistic to 
do conduct tests for atmospheric contaminants at room 
temperatures. 

VIII. CONCLUSION 

In conclusion, we developed a simple thermodynamic 
model to describe the heat transfer mechanisms and gen- 



eration of acoustic waves in photoacoustic Raman spec- 
troscopy for detection of trace amounts of small solid 
particles suspended in air. As it has been pointed out 
since its first application '6'| PARS has several advan- 
tages over other stimulated Raman techniques: it is not 
plagued by a non-resonant background and there is no 
necessity for phase-matching, two very crucial issues in 
Coherent Anti-Stokes Raman Spectroscopy (CARS) for 
example, also PARS directly measures the deposited Ra- 
man energy instead of a fractional energy change. The 
noise in PARS arises primarily from thermal fluctuations 
which may be reduced by controlling macroscopic pa- 
rameters, unlike in other Raman process where the noise 
arises from microscopic quantum fluctuations harder to 
control. Even though thermal noise in a gas is well docu- 
mented, we applied Langevin methods to calculate it for 
our particular physical situation; a reason for doing that 
is that it can be immediately applied to include other 
sources of noises in more detailed studies in the future. 

We applied our model to the specific case of anthrax 
spores and showed that as few as 100 spores in a volume 
of one cubic meter could be detected if thermal noise is 
the main limitation in the system. This promises poten- 
tial development of a PARS-based method to detect such 
harmful contaminants in real time with relatively simple 
means at room temperatures, something that is certainly 
of very strong practical and immediate interest. 
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